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Abstract 

Spontaneous symmetry breaking usually occurs due to the tachyonic (spin- 
odal) instability of a scalar field near the top of its effective potential at = 0. 
Naively, one might expect the field (p to fall from the top of the effective po- 
tential and then experience a long stage of oscillations with amplitude 0{v) 
near the minimum of the effective potential at (j) = v until it gives its en- 
ergy to particles produced during these oscillations. However, it was recently 
found that the tachyonic instability rapidly converts most of the potential 
energy V{0) into the energy of colliding classical waves of the scalar field. 
This conversion, which was called "tachyonic preheating," is so efficient that 
symmetry breaking typically completes within a single oscillation of the field 
distribution as it rolls towards the minimum of its effective potential [Q]. In 
this paper we give a detailed description of tachyonic preheating and show 
that the dynamics of this process crucially depend on the shape of the effective 
potential near its maximum. In the simplest models where V{(l)) ^ —m'^cjp' /2 
near the maximum, the process occurs solely due to the tachyonic instability, 
whereas in the theories —Ac;/)" with n> 2 one encounters a combination of the 
effects of tunneling, tachyonic instability and bubble wall collisions. 
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I. INTRODUCTION 



Since the beginning of the 1970's, spontaneous symmetry breaking (SSB) has been a 
basic feature of all realistic theories of elementary particles. It is discussed in every book on 
quantum field theory, so one might expect the theory of this effect to be well understood. 
However, until very recently this was not the case. 

The standard picture of SSB that many people have in mind looks as follows. The field 
initially stays at the top of the effective potential V{(f)) at = like a ball at the top of a 
hill. Then some small external force pushes it to the right or the left. Even if this force is 
infinitesimally small, it is enough for the field to start falling down in the direction in which 
it was pushed. The field then oscillates near the minimum of its effective potential ^(0) 
at 1 01 = V, and eventually the whole universe becomes filled by a homogeneous scalar field 
\(f)\ = V. As an example of this process one can imagine a piece of ferromagnetic material 
being magnetized under the influence of a very small external magnetic field. 

That is why many authors who studied SSB assumed that the field (p was initially slightly 
displaced from the maximum of V{(f)). Then they studied classical rolling of the field (p from 
this displaced state and the growth of quantum fluctuations on top of the homogeneous 
classical field. It was generally thought that the stage of oscillations of a homogeneous 
classical field with amplitude |A0| ~ v would last for a very long time until it produced 
elementary particles that drained the energy of the classical oscillations. One can talk 
about spontaneous symmetry breaking only when the field settles down at |0| ~ f and 
starts oscillating near this state with an amplitude much smaller than v. |A0| <^ v. The 
most efficient process that was previously known to convert the energy of a homogeneously 
oscillating scalar field into the energy of elementary particles and make the amplitude of the 
oscillations small was parametric resonance 0], but in most cases studied in the literature 
the homogeneous component of the field makes several dozen oscillations before the process 
completes. 

However, in a recent paper by Felder, Garcia-Bellido, Greene, Kofman, Linde and 
Tkachev it was shown that typically spontaneous symmetry breaking completes within 
a single oscillation of the scalar field. One key observation made in |I| was that nobody 
pushes the field from the top of the effective potential in the early universe, so the usual 
picture of a homogeneously oscillating scalar field is incorrect in application to SSB. In those 
cases when the initial value of the homogeneous component of the field is close to zero, 
quantum fluctuations rather than the classical rolling of the homogeneous field dominate 
the dynamics. 

Usually SSB occurs because of the presence of tachyonic mass terms such as —m?(fP'/2 
in the effective potential, so that V" = —m? < 0. Long wavelength quantum fluctuations 
0fc of the field with momenta k < m grow exponentially, 0^ ~ exp{ty/m'^ — k"^). When 
these fluctuations become large they can be interpreted as classical waves of the scalar field. 
Spontaneous symmetry breaking occurs when the total amplitude of these fluctuations grows 
up to V. Because all modes with k < m are growing, SSB occurs locally on a scale somewhat 
greater than m~^. Later on, this scale gradually increases. Inhomogeneities of the scalar 
field absorb some part of the energy ^(0), which suppresses the amplitude of the scalar field 
oscillations. As a result, the field that appears after SSB is relatively homogeneous on 
the scale somewhat greater than m~^, and the amplitude of its oscillations A0(x, t) about 
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the state |0| ~ w is substantially smaller than v. Thus, contrary to naive expectations, a 
prolonged stage of oscillations of a homogeneous component of the scalar fields during SSB 
usually does not exist. 

As we will show in this paper, one reaches a similar conclusion even if the field initially 
has been slightly displaced from the top of the potential (no SSB). With this initial condition, 
the homogeneous background field decays within few oscillations due to the broad parametric 
resonance enhanced by the tachyonic regime. 

The process of rapid transfer of the energy of the scalar field ^(0) into the energy of 
its inhomogeneous oscillations due to tachyonic instability was called tachyonic preheating 
One should distinguish between the tachyonic preheating and spinodal (tachyonic) 
instability, which occurs at the very beginning of this process. The first stages of the 
process of SSB related to tachyonic (spinodal) instability can be studied by relatively simple 
methods. However, very soon the process becomes exceedingly complicated. When the field 
grows sufficiently large, one should take into account nonlinear effects. Oscillations of the 
field can trigger an explosive process of particle production due to parametric resonance , 
which can be especially efficient in our case because of the tachyonic instability. Particles 
(waves) of the classical field produced in this process begin interacting with each other 
(rescattering) . At this stage even advanced methods based on the Hartree or |0] 
approximations fail to describe the situation correctly. In addition, from the very beginning 
of the process there may be production of topological defects, which cannot be described 
by perturbation theory. One might expect that since this is a nonperturbative phenomenon 
it cannot materially affect the process of SSB. As we will see, however, the production of 
topological defects is not a small correction but an important feature of SSB. Topological 
defects, like other inhomogeneities generated by tachyonic instability, drain the energy of 
the scalar field rolling down to the minimum of the effective potential. By doing so, they 
diminish the amplitude of subsequent oscillations of the scalar field. 

There is an extensive literature describing SSB, spinodal instability and the production 
of topological defects during high temperature phase transitions in cosmology . To study 
these issues one should find how the temperature changes in the early universe and use 
numerical methods to find out how symmetry breaking occurs in an expanding universe 
with a time-dependent temperature. Many interesting results have been obtained in this 
direction, see e.g. However, most of these results were strongly model-dependent because 
the answers crucially depend on the ratios between masses of the particles, their coupling 
constants, the temperature of the universe and the rate of expansion. To avoid this problem, 
in this paper we will concentrate on the simplest possibility when the temperature was zero 
from the very beginning and the field was standing on the top of the effective potential. This 
will allow us to study basic features of the process of spontaneous symmetry breaking in its 
pure form without extra complications related to high temperature effects and cosmological 
evolution.^] 

Even in this regime, the theory of SSB remains extremely complicated since for its inves- 
tigation one should go beyond perturbation theory. Fortunately, during the last few years 
new methods of lattice simulations have been developed. They are based on the observation 



A discussion of SSB and tachyonic preheating in the early universe will be contained in 
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that quantum states of bose fields witli large occupation numbers can be interpreted as clas- 
sical waves and their dynamics can be fully analysed by solving relativistic wave equations 
on a lattice [|1^,^. Similar methods were used in |T^-|T6[ in application to sphaleron ef- 



fects, the formation of disoriented chiral condensates, and the problem of topological defect 
production in the early universe. In our paper, which extends the previous work we will 
further develop these methods. 

Usually the main output of lattice simulations is the calculation of correlation functions, 
Wilson loops, etc. A significant advantage of our methods is that the semi-classical nature 
of the effects under investigation allows us to have a clear visual picture of all the processes 
involved. One can really see the process of spontaneous symmetry breaking , which helps 
enormously in understanding the nature of this effect. That is why this paper is accompanied 
by many figures that show the development of symmetry breaking in various models. 

In addition to the simplest models with V^"(0) = —m? < 0, we will study some models 
where the curvature of the effective potential near = is negative, but it vanishes at = 0. 
This happens in such theories as — or —X(f)^. (Potentials of the type of — appear in 
the simplest SUSY motivated models of hybrid inflation.) As we will see, the development 
of tachyonic instability in such models is accompanied by bubble formation and growth and 
bubble wall collisions. Rather interestingly, in these theories bubble formation occurs via 



tunneling even though there is no potential barrier in V{(p) ||T7,18|. Moreover in the theory 
—X(f)^ this process occurs even though there are no instantons describing bubble formation in 
this theory. To understand this process one should use the stochastic approach to tunneling 
developed in |19|. 

Section ^ describes the basic theory of spontaneous symmetry breaking and tachyonic 
instability, focusing particularly on the simplest example of a negative quadratic potential. 
In this section we also discuss the definition of occupation number used throughout the 
paper to describe the growth of fluctuations. Section |ITl| generalizes the theory to a broader 
class of potentials and to the case where the homogeneous field begins displaced from the 
maximum of the potential. Section |^ presents the results of our numerical simulations for 
the simplest SSB model, a single field with a quadratic tachyonic term (i.e. V ~ — m^0^/2). 
Section [V| compares our results with results obtained from perturbative calculations and 
shows how and when the perturbative calculations break down. Section extends our 
numerical calculations to the case of a complex field with a quadratic tachyonic term. 

The next two sections discuss the somewhat more complicated situation that arises when 
the tachyonic mass is dependent and vanishes at = 0. Section |Vli] discusses quartic po- 
tentials {V ~ — A0^/4) and explains how SSB occurs through tunneling and bubble formation 
in such models. Section |Vlli] discusses cubic potentials where the behavior is in some ways 
intermediate between that of the quadratic and quartic cases. 

In the concluding section we summarize our results and briefly discuss their application 
to various cosmological scenarios, including hybrid inflation, new inflation, brane inflation, 
and the recently proposed ekpyrotic/pyrotechnic universe scenarios. Finally there is an 
appendix that provides details on our lattice calculations and lists the parameters used for 
each of the simulations described in the paper. 
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II. TACHYONIC INSTABILITY AND SPONTANEOUS SYMMETRY 

BREAKING 



The simplest model of spontaneous symmetry breaking is based on the theory with the 
effective potential 

V{^) = ^{<P'-vy^-^<P^ + ^^^ + ^, (1) 

where = Af ^ and A ^ 1. V{(f)) has a maximum at = with curvature V" = V^^ = — 
and a minimum a.t (p = ±v. 

The development of tachyonic instability in this model depends on the initial conditions. 
We will assume that initially the symmetry is completely restored so that the field does 
not have any homogeneous component, i.e. (0) = 0. But then (0) remains zero at all later 
stages and for the investigation of SSB one needs to find the spatial distribution of the field 
(j){x,t). To avoid this complication, many authors assume that there is a small but finite 
initial homogeneous background field 0(t), and even smaller quantum fiuctuations 6(j){x,t) 
that grow on top of it. This approximation may provide some interesting information, but 
quite often it is inadequate. In particular, it does not describe the creation of topological 
defects, which, as we will see, is not a small nonperturbative correction but an important 
part of the problem. 

Let us consider equation for the scalar field fiuctuations in the model (0): 

0fc + + V")<Pu = . (2) 

For definiteness, we suppose that the mode functions describing quantum fiuctuations in 
the symmetric phase = at the moment close to t = are the same as for a massless 
field, 0fc = "^^^'^^^^^^^ ■ Then at t = we 'turn on' the term —rn?(j)'^/2 corresponding to 

the negative mass squared —m?. The modes with k = \k\ < m grow exponentially. Initial 
dispersion of all growing fiuctuations with k < m was given by 







W^>-/^-g. (3) 



and the average initial amplitude of all fiuctuations with k < m was given by 

H=^- (4) 
Ztx 



The dispersion of the growing modes at t > is given by 

dk^ ^2tV7^ e2-*(2mt - 1) + 1 





This means that the average amplitude 6(j){k) of quantum fiuctuations with momenta ~ k 
initially was 6(f>{k) ~ k/2n, and then it started growing as e*^™^~'^^. 

To get a qualitative understanding of the process of spontaneous symmetry breaking, 
instead of many growing waves with momenta < m in (13) let us consider first a single 
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sinusoidal wave 50 = A(t) coskx with k m and with initial amplitude A(t) ~ in one- 
dimensional space (so that the average value of (50)^ corresponds to The amplitude of 
this wave grows exponentially until it becomes 0{v) ~ m/\fX. This leads to the division of 
the universe into domains of size 0{m~^) in which the field changes from 0{v) to 0{—v). 
The gradient energy density of domain walls separating areas with positive and negative 
will be ~ k'^dcfP' = 0{m'^/\). This energy is of the same order as the total initial potential 
energy of the field V{Q) = m^/4A. This is one of the reasons why any approximation based 
on perturbation theory and ignoring topological defect production cannot give a correct 
description of the process of spontaneous symmetry breaking. 

Thus a substantial part of the energy V^(0) is transferred to the gradient energy of the 
field when it rolls down to the minimum of V{(f)). Because the initial state contains many 
quantum fluctuations with different phases growing at different rates, the resulting field 
distribution is a gaussian random field with varying spectrum. It cannot coherently give all 
of its gradient energy back and return to its initial state = 0. This is one of the reasons 
why spontaneous symmetry breaking and the main stage of preheating in this model may 
occur within a single oscillation of the field 0. 

Meanwhile if one were to make the usual assumption that initially there exists a small 
homogeneous background field <^ f with an amplitude greater than the amplitude of 
the growing quantum fluctuations 50, so that m/27r <^ < m/ a/A, one would find out that 
when falls to the minimum of the effective potential the gradient energy of the fluctuations 
remains relatively small. In some situations this could lead one to falsely conclude that the 
field will experience many fluctuations before it relaxes near the minimum of V{(j)). To avoid 
this error, we need to perform a complete study of the growth of all tachyonic modes and 
their subsequent interaction without making this simplifying assumption about the existence 
of a homogeneous field 0. 



The tachyonic growth of all fluctuations with k < m continues until y (50^) reaches the 

value ~ v/2, since at ~ v/\/3 the curvature of the effective potential vanishes and instead 
of tachyonic growth one has the usual oscillations of all the modes. Eq. (j^) shows that this 
happens within a time t^, ~ In y , where C ~ 10^. 

A convenient tool for studying this process is the spectrum of the growing quantum 
fluctuations. Rather than using the usual power spectrum |0a;P, however, we find it more 
informative to investigate the occupation number ra^ of produced particles [^. Indeed, in 
situations where the number of particles is well defined (and this always happens at the 
end of the process) the occupation number ra^ is an adiabatic invariant, i.e. it does not 
change during the field oscillations unless some dramatic changes occur to the system. The 
standard definition of the occupation number which was extensively used in the theory of 
preheating 0, and which is valid for > 0, is 

-I 

However, this definition does not work in the tachyonic regime when the effective mass 
squared of the field becomes negative since then Uk = \^k^ + rn^ becomes imaginary. 
Strictly speaking, Uk should not be interpreted as the occupation number of particles during 
the tachyonic regime. One may still formally calculate the function Uk in this regime using 
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either \Jk'^ + \im?\ oi \ k\ instead of uj^ = ^/k? + m? in the expression for Uk whenever rn? < 0. 

The choice between ^J'k^~+~\m?\ and |fc| is arbitrary, and it does not change any of the final 
physical results. The spectra shown in this paper used Uk = \k\ in the tachyonic regime, 
with the exception of Fig. 1 and Fig. 2, where we used Uk = \Jk'^ + \rn?\. The formally 
defined quantity Uk can be interpreted as the occupation number of particles after the end 
of the tachyonic regime, when > 0. Moreover, for all nonvanishing momenta these two 
quantities match very well when one switches from the tachyonic regime to the normal one. 
Indeed, whereas the value of Uk changes during the process, the exponential growth of Uk is 
mainly determined by \(f)k\ and \4>k\, which do not change strongly during the switch between 
the tachyonic regime and the normal one. Therefore the function is very convenient and 
informative during the whole process. When one calculates during the tachyonic regime, 
one can get a good idea of the number of particles that will emerge at the end of this regime 
where the usual particle interpretation becomes possible. In this sense we will interpret the 
function Uk defined above as the occupation number of particles in both regimes. 

An additional caveat of this interpretation is that when one calculates 0^, one does not 
distinguish between the contribution to this quantity from small perturbations and from 
topological defects. As a result, in the presence of topological defects one can somewhat 
overestimate the number of produced particles. The error, however, is not very large, espe- 
cially in theories where instead of domain walls we have strings or monopoles. Moreover, 
eventually topological defects disappear and release their energy in the form of produced 
particles, and the standard interpretation of becomes completely valid. 

The exponential growth of fluctuations during the tachyonic regime can be interpreted 
as the growth of the occupation number of particles with k <^m. Using the estimates given 
above, one can show that Uk for <^ m at the time t^, grows up to 

Hk ~ exp(2mt) = 0(10^) A"^ > 1 . (7) 

The time t^, ~ 2^ In j depends only logarithmically not only on A, but, more generally, on 
the choice of the initial distribution of quantum fluctuations. As we see, for small A the 
fluctuations with k <^ m acquire very large occupation numbers. More importantly, such 
fluctuations will have large amplitude and will be in squeezed state. The general solution 
for these fluctuations will contain two terms, Ae'^^ and Be~^^, but after a short time only 
the growing mode Ae^^ survives. Therefore, independently of the initial phases of quantum 
fluctuations, the only modes with k < m that survive after the beginning of the tachyonic 
regime will coherently grow, and their amplitude will become extremely large. That is 
why these modes can be interpreted as classical waves and can be studied by computer 
simulations using the methods of |10|,|Tl|. 

The dominant contribution to (50^) in Eq. (|^) at the moment t^, is given by the modes 
with wavelength ~ 2nk~^ ~ y/27rm~^ In^/^ (C/A) > where C = O(IO^). As a result, 
at the moment when the fluctuations of the field reach the minimum of the effective 
potential, (Scj)'^) ~ v, the field distribution looks rather homogeneous on a scale I <^U. On 
average, one still has (0) = 0. This implies that the universe becomes divided into domains 
with two different types of spontaneous symmetry breaking, (j) ~ ±f . The typical size of 
each domain is U/2 ~ m^^ In^''^ y, which is slightly greater than m~^. At later stages 
the domains grow in size and percolate (eat each other), and SSB becomes established on a 
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macroscopic scale. 

Of course, these are just simple estimates that should be followed by a detailed quan- 
titative investigation. When the field rolls down to the minimum of its effective potential 
its fluctuations scatter off each other as classical waves due to the interaction. It is 
difficult to study this process analytically, but fortunately one can do it numerically using 
the method of lattice simulations developed in [|l^,|ll|]. 

Before describing the results of our lattice simulations, we would like to discuss the setting 
of the problem, the choice of the initial conditions and some other aspects of tachyonic 
instability in a more general class of theories, including the theories with V{(f)) ~ — 



III. TACHYONIC INSTABILITY IN A MORE GENERAL CLASS OF 
THEORIES AND THE ROLE OF INITIAL DISPLACEMENT 

As we already emphasized, spontaneous symmetry breaking usually occurs due to quan- 
tum fluctuations when the field falls from an exactly symmetric state (f) = 0. However, it 
is still very instructive to find out what happens when the field falls down from some state 
with 00 7^ 0. By doing this, we will get an additional insight in the nature of tachyonic 
instability. We will also be able to compare our results with the results of earlier works on 
spontaneous symmetry breaking. 

Consider the behavior of the fluctuations 6(f) with momentum k <^ m. An important 
observation is that these fluctuations satisfy the same equation of motion as 0: 

0fe = -y"(0)0,. (8) 

A general solution of this equation for V = — m^0^/2 is 0fc(t) = aie™* + 026""^*. Similarly, 
for one has = 6ie™'* + 62e~™*. At t ^ rrT^ only the growing mode survives, and the ratio 
0fc/0 becomes constant. This rule holds for other types of tachyonic potentials as well. Thus 
one can investigate the amplification of the long wavelength perturbations of the scalar field 
in a very easy way. Instead of solving equations for 0^ in a time-dependent background 
0(t), one can find how 0(t) changes in time. We will use this trick here and in the next 
section. 

Consider the theory 

V(0) = - A072 . (9) 

Suppose the field begins rolling down from 0o. Energy conservation implies that 

0V2 - 02/2 = V(0o) - V{<P) . (10) 

We will assume for simplicity that in the beginning, at = 0o, the field moves with the 
same velocity as if it were falling with vanishing total energy from = (this assumption 
does not make any difference for motion at ^ 0o). Then one has 

02/2 = -r(0) = A0"/2 . (11) 

Thus 
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(/. = yA0"/2 . (12) 

The solution is 

>„--VAt(!i^))-*. (13) 

The most important result is Eq. (|T^. In a more general case of nonvanishing total energy 
we have 



= ^A((/." - 0^) , (14) 

where (po is an initial field value where 4> = 0. It implies that the tachyonic fiuctuations with 
small momenta in the long time limit (0 ^ 0o) grow as follows: 

0, = c0"/^ (15) 

where C is some constant. 

This means, in particular, that in the theory with the potential —0^ the long wavelength 
fiuctuations grow just like the field itself, ^ = const. Meanwhile for the theory with —0^ the 

fiuctuations grow faster, ^ ~ 4>^^'^, and for the theory —0^ they grow even faster, ^ ~ 0. 

Returning to the theory —^0^ + ^0^ + we find that the potential is tachyonic 
{V" < 0) for < < v/VS, and it can be approximately represented as —^0^ for 
< ^ t'/2. When the field grows from 0o <^ f to f / -\/3, the speed of the field grows 
from m0o to ^mv. Consequently, the amplitude of density perturbations grows by a factor 
~ and the occupation numbers of particles for A; <^ m grow by a factor O (^)- 

Clearly, one has the largest amplification if one starts as close to 0o = as possible. 
However, if 0o <C where ^ is the average amplitude of the long wavelength quantum 
fiuctuations with momentum k < m (which grow almost as fast as the homogeneous mode), 
then the development of 0o gives no information on the process of spontaneous symmetry 
breaking. In this case instead of 0o one would need to study all growing modes with k < m, 
just as in the case 0o = 0. This is what we are doing in the main part of this paper. 

Equation for fiuctuations in the model V{(j)) = —^0^ + f 0^ + ^ is 

0fc + {k^ -m^ + 3A02) 0fc = , (16) 

This equation should be solved simultaneously with the equation for the background field 
0(t) 

- m20 + A03 = , (17) 

Equation (|l^) is the Lame equation [^]. Its solutions depend on the dimensionless parame- 
ters \/X(f)o/m. In the context of the chaotic infiationary model, where the field 0(t) is rolling 
from its large initial value ~ Mp, this parameter usually was taken to be large ||2^. In the 
context of the theory of spontaneous symmetry breaking we are dealing with the opposite 
case when 0o is close to zero. 
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The description of the growth of perturbations in the model —^0^ + ^4''^ + ^ is a 
straightforward generahzation of the theory of parametric resonance in the model which 



has been studied using the stability/instability chart of the Lame equation |]20[. The presence 
of the negative mass term adds an additional instability band at k ^m. The characteristic 
exponent /i in this new zone is significantly greater that in the higher zones because of the 
tachyonic effect. Thus, the tachyonic parametric resonance will be dominant. 

When the field rolls towards the minimum of V{(f)), the occupation numbers Uk, calculated 
from the solutions 0^ of equation (|l6l), become large. However, for (po > ^ the field 
fluctuations do not grow large enough to dominate the energy density immediately after the 
rolling to the minimum of the effective potential. To find out what happens in this case, we 
will describe, as an example, the evolution of the occupation numbers of the modes (pk with 
different momenta k in the model (|ID with A = 10~^ if the field rolls from 0o = O.Olv (which 
is larger than ^ in this model). 

Consider first a mode 0^ with k ^m. In the beginning, when the field rolls from (p = (pQ 
to (p = f/v^, this mode grows faster than any other fluctuations, just as we expected, see 
Fig. ^ During this time interval, the occupation number becomes ~ e^, which is in good 
agreement with our estimate O ~ 2 x 10^. Then the field reaches the bottom of the 

effective potential, goes somewhat beyond this point, bounces back, and again approaches 
the tachyonic region (p < Until the field becomes smaller than v/\^, the occupation 

number of particles with k <^ m does not change much. But then it decreases almost to 
the same value from which we started our calculations. What happens is that the solution 
for the fluctuations has two modes, the growing one and the decaying one. When the field 
bounces, fluctuations either grow or decay depending on the phase with which they re-enter 
the tachyonic regime. Therefore even though the modes with -C m grow fast on the way 
down, they also decay fast on the way up, as shown in Fig. |l]. 

Meanwhile for the modes in a rather broad interval of k, from k ~ 0.3m to ~ 0.6m, 
the modes continue their growth when the field oscillates. Figure ^ shows the growth of 
during three consecutive oscillations of the field (p. As we see, during each full oscillation 
the occupation numbers grow e^^ times. Thus during n oscillations the occupation numbers 
should grow e^^" times. 

This is an incredibly fast growth. It occurs much faster than the usual parametric 
resonance in theories with m^ > 0. Clearly, this process can rapidly convert all the 
energy of the homogeneous field into the energy of classical colliding waves, and at this 
stage the only reliable way to study the process is to use numerical simulations. 
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FIG. 1. Evolution of the occupation numbers for the fluctuations with k m in the model 
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FIG. 2. Same as in Fig 1 for k = 0.5m. 



In the theory — m^0^ the fluctuations grow as fast as the scalar field, so if one begins 
with a homogeneous field with ^ -C 1, then on the way down to ~ w this field distribution 
remains relatively homogeneous. However, when the field rolls back towards (p = (po, the 
inhomogeneities with k ~ 0.5m continue growing. For example, if one takes 0o ~ 10~^w, 
quantum fluctuations (which have initial amplitude only one order of magnitude smaller than 
00 in this model) grow almost 10^ times when the field falls down from = 0o and returns 
back. The amplitude of inhomogeneities after the return becomes approximately three 
orders of magnitude larger than 0o, which means that the homogeneity becomes completely 



11 



destroyed. At this stage (and in fact much earher) one can no longer study the evolution of 
quantum fluctuations as if they were small deviations on a homogeneous background. When 
the field falls down to the minimum of the effective potential again, it becomes divided into 
large colliding waves. One cannot study the evolution of such a system using perturbation 
theory. 

In the theories V{(j)) ~ —(f)"^ with n > 2 the situation may be even more interesting and 
the growth of the occupation number Uk for small k occurs even faster. For example, in the 
theory — long wavelength fluctuations grow as 0^ (and the occupation numbers grow as 
0^). Therefore when the field grows from 0o to v, the ratio ^ grows by a factor of O (^-^^ . 
This means that the field may become very inhomogeneous on its way down even if initially 
it was very homogeneous. 

The average initial amplitude of tachyonic fluctuations in the theory —X(j)'^/4 at 0o 7^ 

is given by 6(j) ~ "^j^ ^ ~ Initial level of inhomogeneities was given by ^ ~ -C 

1. When the field reaches some value v ^ (f)Q, the ratio ^ grows and becomes ^ ~ 



Thus, if the rolling of the field begins at a very small value of the field 0O; or if it 

continues long enough, so that ^ becomes greater than 1 / V^, the field becomes completely 
inhomogeneous on its way down. 

Moreover, if one considers a theory such as e.g. ^(0) = — A0^/4 + A0^/f^, which has 
a minimum at = f , then in such theories, just as in the theory — m^0^, there are some 
modes with k ~ V^0o whose amplitude grows both on the way down and on the way 
up. For these modes the degree of inhomogeneity rapidly grow with each oscillation. The 

occupation numbers grow approximately as (^^^ during each full oscillation, so that after n 

oscillations the occupation numbers of the particles with momenta ~ "\/A0o become as large 

as . That is why it takes only one or two oscillations before the oscillating scalar field 

becomes inhomogeneous and the first stage of preheating related to the tachyonic instability 
completes. 

The simple rules derived above explain the extraordinary efficiency of tachyonic preheat- 
ing. However, one can apply these rules only at the beginning of the process, when one can 
neglect the backreaction of created particles. That is why we needed to perform computer 
simulations which took the effects of backreaction into account. 



IV. LATTICE SIMULATIONS OF SPONTANEOUS SYMMETRY BREAKING IN 
THE THEORIES WITH V{(p) = -M'^(t?/2 + A0^/4 



A description of our method |[T3|T3I application to this problem is given in the ap- 
pendix. 

There are several complementary ways one can represent the results of our calculations. 
One of the best ways to do it is to study the probability distribution function P(0, t), which 
is the fraction of the volume containing the field at a time t. At t = we begin with the 
probability distribution concentrated near = 0, with the quantum mechanical dispersion 
(H), and then we follow its evolution; see Fig. R|. 
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FIG. 3. The process of symmetry breaking in the model (|l|) for A = 10~^. The values of the 
field are shown in units of v, time is shown in units m~^. For each moment of time, we also show 
the occupation numbers nt (the lower part of each panel), with k measured in units of m. At t = 
one has Uk = 0, as in the usual quantum field theory vacuum. In the beginning of the process the 
occupation numbers nt grow exponentially for A; < m (A; < 1 in the figure), but then this growth 
spreads to k > m because of domain wall formation and collisions of classical waves of the field 
(p. Within a single oscillation the occupation numbers for k <^ m grow up to ~ 10^, which is in 
complete agreement with our estimate nk ~ 10^A~^, Eq. (0). The spectrum rapidly stabilizes, but 
it is not thermal yet, and the occupation numbers remain extremely large. Thermalization takes 
much more time than spontaneous symmetry breaking. 

In the beginning quantum fluctuations are very small, and the probability distribution 
P(0, t) is very narrowly focused near = 0. Then it spreads out and shows two maxima 
that oscillate about (p = ±v with an amplitude much smaller than v. 

As we see from Fig. ^, the two maxima never come close to the initial point = 0, 
which implies that symmetry becomes broken within a single oscillation of the distribution 
of the field 0. To demonstrate that this is not a strong coupling effect, we show the results 
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for the model (|l]) with A = 10~^. We obtained similar results for A = 10~^. Note that 
only when the distribution stabilizes and the domains become large can one use the stan- 
dard language of perturbation theory describing scalar particles as excitations on a (locally) 
homogeneous background. That is why the use of the nonperturbative approach based on 
lattice simulations was so important for our investigation. 

One may wonder why the distribution is slightly asymmetric, and why after symmetry 
breaking there are still many points at -C v. The answer is that after SSB space becomes 
divided into domains with cf) ~ ±f . Domains are large, and their size gradually grows after 
SSB because large domains "eat" the small ones. Eventually in any finite size box there 
will remain just one domain, i.e. the distribution will become completely asymmetric. The 
points with |0| ^ f correspond to domain walls. 

In this series of simulations we made a cut-off in the spectrum of initial fluctuations at 
k > m. The reason is that only the modes with k < m from the very beginning experience 
exponential growth and behave as classical fields. We checked, however, that the results of 
the simulations remain qualitatively the same if one makes a cut-off at k ^ m. 

The process of thermalization takes much longer than spontaneous symmetry breaking 
|21| . Indeed, the standard thermal distribution is given by the well known equation Uk = 
i^e^k/T _ i)-i. At the moment when all the energy 1^(0) = ^ is transferred to the thermal 
energy ~ T^, the temperature rises up to T ~ mA~^/^, and the occupation numbers aXk 
become ~ (e™/-^ — 1)^^ ~ T/m ~ A^^/^. In particular, for A ~ 10~^ one would have 
nk<m = C'(IO), which is 5 orders of magnitude smaller than the results of our calculations. 

Thus, the occupation numbers should drop down dramatically before full thermalization 
is achieved. This may happen only if the total number of particles becomes many orders of 
magnitude smaller (particle cannibalism). If one considers only those interactions that pre- 
serve the total number of particles (scattering 2 — 2), one may achieve a kind of temporary 
thermal equilibrium with a nonvanishing effective chemical potential of particles 0. This 
is what we see in our calculations when the occupation numbers gradually approach some 
quasi-equilibrium asymptotic limit. Since there is no real particle conservation in this the- 
ory, eventually the effective chemical potential will vanish, and the true thermal equilibrium 
with rifc = (^e'^k/T _ -^^-i ^jjj reached. But this process takes much greater time than the 
time required for spontaneous symmetry breaking. 

To provide a visual picture of the distribution of the scalar field, we show the growth 
of fluctuations in a two-dimensional slice of 3D space in this model in Fig. ^. Maxima 
correspond to domains with > 0; minima correspond to domains with < 0. The third 
image corresponds to the first one half of an oscillation, just like the third panel in Fig. ^. As 
we see the universe at that moment is already divided into domains with ~ ±v. The initial 
size of each domain is somewhat greater than m~^. Inside each domain the deviation from 
= 1^1 is much smaller than v. This confirms our conclusion that spontaneous symmetry 
breaking occurs within a single oscillation. 
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The original domain structure can change within a time O(10m~^) because of domain 
wall collisions and domain expansion. Gradually, the size of each domain grows and the 
domain wall structure becomes more and more stable, as we see in the last two images of 
Fig. |. 

If one continues the calculation for a much longer time, one can see much more clearly 
the formation and growth of domains with (j) = ±v. In the beginning these domains are 
small, but then they "eat" each other and grow. To illustrate this process we performed 
simulations in a 2D box of size 1024 x 1024. This allowed us to perform the calculations 
out to a much greater time and see domain formation on a much greater scale; see Fig. |^. 
One should note that the process of rescattering of particles produced during preheating in 
2D is somewhat different from that in 3D. However, the tachyonic instability is the same in 
both cases and the process of domain growth is qualitatively similar. 



t = t = 100 




FIG. 5. Formation of domains in the process of symmetry breaking in the model (Q). 
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V. COMPARISON WITH THE USUAL PERTURBATIVE APPROACH 



Before going any further, let us discuss the difference between our methods and the usual 
approach based on perturbation theory. In the usual approach one solves a self-consistent 
system of equations for the homogeneous scalar field and for the variance (50^), where 
(50^) is generated due to particle production by the oscillating field 0, see e.g. 

In this approach one would expect that the effective mass of the field is given by 

= — + 3A((50^). This is the standard approach used, in particular, in the theory of 
the high temperature cosmological phase transitions 0. 

Let us see what would happen if we naively applied this method to our problem. Our 
investigation shows that within a single oscillation the variance of the scalar field grows to 
(50^) ^ f^; see Fig. ^. A more detailed investigation shows that soon after the beginning of 
the process the value of (50^) is not much different from |f ^. This would seem to imply that 
immediately after spontaneous symmetry breaking the symmetry becomes restored again, 
because the effective mass squared of the field becomes positive: m| ~ — m^ + |Af ^ ~ |m^. 

However, our numerical calculations clearly demonstrate that this is not the case. So 
what is wrong with the standard perturbative approach? 

The point is that in order to investigate the tachyonic instability one should study the 
local field distribution on a scale comparable to instead of the field distribution averaged 
over the whole universe. In our scenario the universe becomes divided into different domains 
of size greater than 0(m"^). If one wants to study the process of spontaneous symmetry 
breaking, then instead of finding the average values (0) and (50^) over the whole universe 
one should find the average value of inside each domain. After that, one should calculate 
the variance (50^), where 50 is the /oca/ deviation of the field from its average value inside 
each domain. But even this will give only partial information about the process. That is 
why in addition to finding the probability distribution and the occupation numbers (Fig. 

we have shown the spatial distribution of the field (Fig. |^). 

This need for local averaging is an important issue that was overlooked in many recent 
works on preheating, as well as in some works on the backreaction of long wavelength infla- 
tionary quantum fluctuations on the speed of expansion and the average energy-momentum 
tensor of matter. We are not saying that the calculation of averages such as (0) and (50^) 
over the whole universe is not useful. For example, it is quite informative in the theory 
of high temperature phase transitions, where the typical contribution to (50^) occurs due 
to the short wavelength fluctuations with the wavelength T~^, which is much smaller than 
m^^ at the time of the phase transition 0. However, one should be extremely careful using 
averages like (50^) over the whole universe in situations where a substantial contribution 
to (50^) is given by fluctuations whose wavelength is greater than the typical length scale 
of the problem. It does not matter how accurately one calculates such averages, whether 
one works in the Hartree approximation or in the approximation, as in In the 
case described above we calculated (50^) very accurately using our lattice simulations. This 
method takes into account not only the usual backreaction effects that could be studied in 
the Hartree or approximations, but also effects of rescattering of produced particles. 
Still we have seen that a naive use of our results would lead to an incorrect conclusion that 
symmetry becomes restored immediately after it breaks down. 

If one has to use perturbation theory in situations when the infrared contribution to 
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b"^) is substantial, the occupation numbers are large and the results allow a semi-classical 
interpretation, one can avoid the problem discussed above if one rearranges perturbation 
theory in a nontrivial way. For example, if one studies effects on a length scale /, one may 
consider all fluctuations on larger scales as a nearly homogeneous classical field background 
and ignore the contribution of these fluctuations to (50^). In the context of the theory of 
preheating, this issue was discussed in Section X of Ref. 0]. In inflationary cosmology a 
similar approximation constitutes the basis of the stochastic approach to inflation P^p5| ,|5|. 
An alternative method is to use numerical simulations that can bring us more detailed 
information about the process. This is the method we use in our paper. 

To compare our method with the more traditional perturbative approach assuming initial 
displacement of the field 0o 7^ 0, we performed a series of simulations for different values of 
00 exceeding the level of the long wavelength quantum fluctuations ~ For 0o ~ 10~^w 
the results did not differ much from the results that we obtained for 0o = in the previous 
section. The distribution of the field never returned back to the vicinity of = 0, and the 
process of spontaneous symmetry breaking occurred within a single oscillation. Note, that 
our calculations were performed for A ~ 10~^, so that 0o ~ 10~^f ^ ™ - ■ 



2-K 2-K ■ 

For a much greater initial displacement the results were somewhat different, but still for 
00 <^ f we have found that the regime of homogeneous oscillations completely disappeared 
after a couple of oscillations. Consider, for example, the case 0o = O.lu; see Fig. |[ As 
one might expect, in this case the final probability distribution is entirely concentrated at 
~ +f , and one can check that no topological defects are produced at the end of the 
process. But even in this case we have found that the process completes very fast, after the 
second oscillation. 

This result differs from the results of investigation of the same model in 0], where it was 
claimed that the oscillations of the homogeneous component of the field in this model con- 
tinue for a long time with the amplitude comparable to f , even if one starts with 0o -C 
The reason for the disagreement is very simple. First of all, in the investigation of sponta- 
neous symmetry breaking in the theory —^0^ + ^0^ + ^ in the initial displacement 
of the field 0o was chosen two orders of magnitude smaller than the level of the long wave- 
length quantum fluctuations with < m, 50 ~ In this case SSB appears not because 
of the growth of the homogeneous component of the field, but because of the generation of 
fluctuations with k < m. In such a situation investigation of the homogeneous component 
of the field does not give much information about spontaneous symmetry breaking. 

Moreover, as soon as the combined amplitude of all fluctuations with k < m (i.e. 



'^)k<m) becomes comparable to f , which happens much earlier than the homogeneous 
component of the field reaches the minimum of the effective potential, the universe becomes 
divided into domains with colliding walls. At this moment the standard perturbative ap- 
proach completely breaks down. It does not describe rescattering of produced particles, 
collisions of classical waves of the scalar field, and dynamics of topological defects. In this 
regime equations describing the evolution of the homogeneous component of the field de- 
rived in the Hartree approximation (or in 1/N approximation) become inapplicable. That is 
why in our work we studied SSB using a combination of analytical investigation and lattice 
simulations. 
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FIG. 6. The process of symmetry breaking in the model (||) for A = 10^^ in the case where 
the field (/) was initially displaced from = by (/>o = 10^^?;. As we see, in this case spontaneous 
symmetry breaking takes two oscillations to occur, and the occupation numbers are smaller than 
in the case of falling from </> = 0. 
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VI. SPONTANEOUS SYMMETRY BREAKING IN THE THEORY OF A 

COMPLEX FIELD 



In the previous section we studied symmetry breaking in a theory (|TD describing a one- 
component real field (f). One can perform a similar investigation for the theory of a multi- 
component scalar field 0i with the potential (|l|), simply replacing 0^ with Figure 
1^ illustrates the dynamics of symmetry breaking in the model (Q) with a two-component 
scalar field = (0i -|- i02)/v^- It shows the probability distribution P(0j,t), which is the 
fraction of the volume containing the field at a time t. 



t=0 t=9 




FIG. 7. The process of symmetry breaking in the model ([l|) for a complex field 
(f) = -^{4>i + *02)- The field distribution falls down to the minimum of the effective potential 
at 1 01 = V and experiences only small oscillations with rapidly decreasing amplitude |A0| <C v. 

In the beginning, the probability distribution is concentrated near = 0, with the 
quantum mechanical dispersion (^. Then the probability distribution spreads out, and 
after a single oscillation it stabilizes at |0| ~ which corresponds to SSB. The standard 
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approximation representing the scalar field as a homogeneous background field with small 
fiuctuations does not work at any stage of the process. 

A detailed investigation of the spatial distribution of the field shows that after the 
first oscillation the scalar field can be represented as a collection of classical waves oscillating 
near \(f)\ ~ v with an amplitude smaller than v/2. Thus SSB indeed occurs within a single 
oscillation of the field distribution. A small but nonvanishing height of the histogram in Fig. 
1^ at = is due to the presence of strings that have = at their cores. 

Fig. P shows the occupation numbers of produced particles. During the first oscillation 
these numbers grow up to 10^ - 10^ for k < m [k < 1 in the figure). Then the occupation 
numbers a.t k < m slightly decrease, whereas the occupation numbers aX k > m begin to 
grow. Complete thermalization takes a very long time. 




FIG. 8. Occupation numbers of particles produced during tachyonic preheating in the 
model of a complex scalar field (j) with effective potential V = —m'^(j)*(l) + X{(p*(j))'^ with A = 10~^. 
In the beginning (lower curves), Uk grows for /c ^ m (/c ^ 1 in this figure), but then eventually this 
growth spreads to larger k. 

In the model of a complex scalar field, instead of domain walls one has strings that are 
produced when the field falls down; see Fig. |^. The whole process of string formation occurs 
within a single oscillation. After that the new long strings are not formed. Sometimes small 
string loops appear and disappear because of occasional large fluctuations of the scalar field. 
If symmetry were broken and then restored again when the field distribution moves back 
to = 0, we would see strings being "melted," and then a completely new set of strings 
would appear. Meanwhile, our simulations show that the large scale string distribution is 
formed as the field first rolls down to the minimum of the effective potential. During the 
subsequent oscillations the strings formed in the beginning of the process do not disappear 
and are not replaced by new ones; instead they experience only gradual evolution. This 
confirms our conclusion that symmetry breaking is achieved within a single oscillation. 

Just as in the case of the one-component scalar field, perturbative methods of investiga- 



21 



tion of this theory cannot describe formation of topological defects and scattering of classical 
waves produced by the tachyonic instability. Therefore such methods break down within 
the first oscillation of the field distribution. 

t= 11 



100 




FIG. 9. Strings produced after one half of an oscillation in the model (Rl) for a complex field 0. 



VII. QUARTIC POTENTIAL 



The process of SSB will occur in a somewhat different way in theories where the curvature 
of the effective potential near its maximum depends on 0. For example, one may consider 
the Coleman- Weinberg model, which was the basis for the first version of the new inflation 
scenario 
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(18) 



This potential has a maximum at = and a minimum ai ip = v. At small the effective 
potential looks hke — j- with an effective coupling constant A(0) = Alog^. 
Another interesting example is the toy model 



A 



6 ^4 



V = 0^ + + 

r 6 12 



(19) 



An important feature of such potentials is that the tachyonic mass m^(0) = \^"(0) 
vanishes at = 0. Therefore the simple arguments based on the tachyonic growth of small 
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quantum fluctuations do not apply here. The decay of the symmetric phase in such models 
occurs via tunneling and the formation of bubbles. Historically, this was the flrst example 
of a theory where tunneling occurs between two states (0 = and 7^ 0) even though there 



is no barrier separating these states [0; see also ||T8 



To study symmetry breaking in these models one should flrst consider the growth of the 
fleld d) in the model 



(20) 



and then see what happens when one adds extra terms that stabilize the potential. 

The tunneling trajectories (instantons) with minimal action possess the 0(4) symmetry 



of Euclidean space p5|. The Euclidean equation for 0(4) symmetric tunneling is 



+ 



^'(0) • 



with the boundary conditions 0(r = oc) = v and 0'(O) = 0. Here ^'(r) 
Xi are the Euchdean coordinates, i = 1,2,3,4. 

Equation (^) in the theory ( pOD has a family of solutions p6| , P!7 



dr ' 



(21) 



xf; the 



0(r) 




P 



(22) 



where p is arbitrary. Note that the value of the scalar fleld in the center of the bubble 
depends on p: 



0(0) 



2^2 

TTp 



(23) 



The corresponding Euclidean action does not depend on p. 



27r^ / rM^(0')' + V^(0))c^r 



(24) 



The probability of bubble formation per unit four- volume can be estimated by the expression 



P p ^ exp 



3A 



A>^(0)exp 



3A 



(25) 



The probability of tunneling in the Coleman- Weinberg theory (|T8|) can be estimated by 
this equation if instead of A one uses the effective coupling constant Alog^. Tunneling is 

2 

not strongly suppressed at Alog^ ~ 1. This means that tunneling occurs to a point with 

exponentially small (p: ~ f e~'^^^, with C = 0(1). 

On the other hand, in the model (|T^) the effective coupling constant A and the factor 
exp (— suppressing the tunneling do not depend on 0, whereas the subexponential factor 
in the expression for the tunneling probability (^) is greater for large 0. Thus in this model 
tunneling may occur to relatively large 0. 
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The bubbles that appear after the tunnehng are described by eq. ( p2D if one understands 
by its Minkowski counterpart — t^: 




t2 + p2 



(26) 



Such bubbles have symmetry 0(3, 1). When the bubble appears (at t = 0), the field takes 
its maximal value 0o at the center of the bubble, 0o = f \/f- Then it grows, and becomes 
infinitely large at t = p = = Here (0) = V"{<P) = 

Of course, in realistic models like ([T8|) and (^) the field does not grow indefinitely large. 
It reaches the minimum of the effective potential at = ±v and then it begins oscillating 
there. Meanwhile quantum fluctuations may grow on top of the smooth instanton solution. 
The investigation of these oscillations and bubble wall collisions is a complicated problem 
that can be studied numerically. Fortunately, the behavior of the oscillating fleld prior to 
the bubble wall collisions and neglecting quantum fluctuations can be studied analytically 
by making a certain change of variables. 

Indeed, it is known that in properly chosen coordinates the interior of each bubble looks 



like an open universe filled by a homogeneous scalar field . One can show that during 
the main part of the first oscillation of the field the radius of curvature (scale factor) of this 
open universe is 0{p), which leads to expansion of the open universe with Hubble constant 
H ~ ~ m0(0o).0 This introduces the damping term 3H(p to the equation of motion 
of the scalar field, which gradually diminishes the amplitude of its oscillations. Suppose 
that the tunneling occurs to 0o <^ "W, as in the Coleman- Weinberg model. Then during the 
main part of the first oscillation the effective mass of the field remains much greater than 
H ~ m<^(0o); so in the limit ^ — >■ one can neglect the effect of expansion of the open 
universe on the amplitude of the oscillations. Later on, the Hubble constant in the open 
universe bubble becomes even smaller and its damping effect on the oscillations becomes 
even less significant. As a result, the amplitude of oscillations of a homogeneous scalar field 
for a long time remains almost unchanged. 

From the point of view of an outside observer using the usual coordinates x and t this 
means that the field in the center of the bubble oscillates for a long time with amplitude 
0{v), sending spherical waves of the same amplitude in all directions. The amplitude of each 
wave is a function of — t^, which means that they propagate with a speed asymptotically 
approaching the speed of light and their amplitude does not depend on their distance from 
the center of the bubble.0 



^In this paper we are neglecting the overall expansion of the universe caused by the energy density 
of the scalar field. This effect will be considered in a separate publication. Here we consider 
"expansion" as it is seen inside the bubble in the coordinate system in which the interior of the 
bubble looks like a homogeneous open universe. This is not a physical expansion but a consequence 
of the choice of the coordinate system in flat space where it is more convenient to study the bubble 
motion. 

^One could speculate about the possibility of sending a signal to aliens using such waves with 
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Thus instead of the naive picture of a bubble consisting of a single spherically symmetric 
shell (which would be a correct picture in the thin- wall approximation), one has a series of 
waves of almost equal amplitude following each other; see Fig. Reheating in this model 
occurs due to a combination of different effects. First of all, particles are produced during 
the collision of waves produced by different tunneling events. But they are also produced 
due to the tachyonic instability, as well as by the oscillations of the scalar field inside each 
bubble. 




FIG. 10. Field values on a partial 2D slice through the lattice in the model 

V = J log ^ ~ ^ + • "^^^ process of symmetry breaking occurs due to tunneling and 
bubble formation. After the tunneling, the bubble grows, and the field inside it begins to oscillate. 
If the tunneling occurs from (/> = to 00 *C w, the amplitude of oscillations remains large for a 
long time, and instead of the usual picture of a single bubble wall propagating in all directions 
one has a series of propagating waves with amplitude comparable to v. The figure shows a half of 
such bubble, which appears after the tunneling to (po = .02v. Cutting the bubble in half allows us 
to see that the amplitude of oscillations decreases rather slowly, just as we expected. Because the 
bubbles in this model take an exponentially long time to form, in our simulations we did not start 
at (/> = and wait for one to appear, but rather started using the analytic form of the instanton as 
our initial conditions. 



an amplitude that does not decrease with the distance. The problem is that these waves are only 
possible as a result of vacuum decay, which first kills those who send the signal and then those 
who receive it. 
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We should make some comments here. First of all, if the tunneling occurs to very 
small values of 0, quantum fluctuations produced due to tachyonic instability inside the 
0(3, 1) symmetric bubble may completely distort the shape of the bubble during the field 
oscillations. Within few oscillations, tachyonic preheating creates colliding waves inside the 



bubble; see Fig. 11 



Note that each tunneling event produces an exponentially large sphere filled either by a 
positive field oscillating around (p = v with a slowly decreasing amplitude, or by a negative 
field 0, oscillating around = —v. In both cases SSB occurs within a single oscillation 
within each bubble, and then finally the field relaxes near ±v due to a combined effect of 
the amplitude decrease because of the bubble expansion, and the development of tachyonic 



preheating, as in Fig. |Tl 



Finally, we should say that if the tunneling occurs to extremely small values of 0, or if 



it does not occur for a long time, one may obtain an inflationary regime pj]. Tachyonic 
preheating in this regime will be discussed in a subsequent publication M. 





FIG. 11. The process of symmetry breaking in the model V = j (^0^ log ^ ~ ^ + \j taking 
into account quantum fluctuations in the instanton background for A ~ 10"^. As we see, quantum 
fluctuations lead to a growing asymmetry and decoherence of the oscillations due to the tachyonic 
preheating inside the bubble. Thus preheating in this model occurs due to a combined effect of 
bubble wall collision and tachyonic preheating. The latter mechanism is especially efficient if the 
tunneling occurs to 0o ^ v. 



26 



VIII. CUBIC POTENTIAL AND STOCHASTIC APPROACH TO TUNNELING 



Another important example of tachyonic preheating is provided by the theory 

V = -^v<p' + ^<p^ + . (27) 

This potential is a prototype of the potential that appears in descriptions of symmetry 
breaking in F-term hybrid inflation ||2^j2^ . 

The development of instability in this theory presents us with a new challenge. The 
curvature of the effective potential at = in this theory vanishes, which means that, 
unlike in the theory —m?(fP' (|I|), infinitesimally small perturbations in this theory do not 
grow. On the other hand, unlike in the theory —Xcf)'^ ([20|), there are no instantons in this 
theory that would describe tunneling from = 0. Thus, in the theory — Aw0^, which occupies 
an intermediate position between — m^0^ and — A^"^, both mechanisms that could lead to 
the development of instability do not work. Does this mean that the state = in this 
theory is, in fact, stable? 

The answer to this question is no; the state = in the theory — At'0^ is unstable. 
Indeed, even though (0) initially is zero, long wavelength fluctuations of the field are 
present, and they may play the same role as the homogeneous field in triggering the 
instability. 

Eq. (^) implies that scalar field fluctuations with momenta k <^ have initial amplitude 
) ~ o-%. Thus the short wavelength fluctuations with momenta k > ko live on top of a 



long wavelength field with an average amphtude 50rms(^o) ~ y (^0^) ~ ivfc' 

The curvature of the effective potential V" = \mlfi\ at ~ 50rms(fco) in the theory 
is given by — 2Av50rms(^o) ~ Consider fluctuations with momentum k somewhat 

greater than /cq, so that the amplitude of the long wavelength field 50 does not change 
significantly on a scale k~^. Short wavelength fluctuations with k = Cko with C somewhat 
greater than 1 will grow on top of the field ~ 50rms(fco) if k"^ ^ l^^efrl ~ 

Taking for definiteness C \/2, one may argue that fluctuations with k ^ ^ may enter 
a self-sustained regime of tachyonic growth. Small fluctuations rapidly grow large, which 
justifies using semi-classical methods for the description of this process. The average initial 
amplitude of the growing tachyonic fluctuations with momenta smaller than ^ is 

50rms ~ (28) 

These fluctuations grow until the amplitude of 50 becomes comparable to 2v/3, and the 
effective tachyonic mass vanishes. At that moment the field can be represented as a collection 



of waves with dispersion y (50^) ~ v, corresponding to coherent states of scalar particles 

with occupation numbers ^fc ~ (^) ^ 1- 

A more accurate investigation shows that the initial value of the field is few times greater 
than 50rms ~ ^ (see below), and therefore the occupation numbers will be somewhat 
smaller, 

Uk ~ 0(10) A-2 . (29) 
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Because of the nonlinear dependence of the tachyonic mass on 0, a detailed description 
of this process is more involved than in the theory (|l|). Indeed, even though the typical 
amplitude of the growing fluctuations is given by (p8D, the speed of the growth of the fluc- 
tuations increases considerably if the initial amplitude is somewhat bigger than (^). Thus 



even though fluctuations with an amplitude a few times greater than ( pSj) are exponentially 
suppressed, they grow faster and may therefore have a greater impact on the process than 
fluctuations with amplitude 

Low probability fluctuations with 5(f) ^ 50rms correspond to peaks of the initial Gaussian 
distribution of the fluctuations of the field 0. The theory of the 3d random Gaussian fields 



is well developed ||3^. Its statistical properties are determined by the spectrum One 
of the most interesting features of the Gaussian field is statistics and the shapes of the high 
peaks of the field distribution. Such peaks tend to be spherically symmetric. As a result, the 
whole process looks not like a uniform growth of all modes, but more like bubble production 
(even though there are no instantons in this model). A simple physical interpretation of 
the inhomogeneous fragmentation of the field is based on the fact that the interaction 
—\v(j)^ corresponds to attraction between the fluctuation modes. As a result, the seed 
inhomogeneities (the peaks of the initial random distribution) will be amplified due to the 
nonlinear interaction of the fluctuations. A well known example of this type of instability is 
gravitational instability of matter in the universe. 

To study the growth of fluctuations in a more detailed way, one may use the stochastic ap- 
proach to tunneling and bubble formation developed in . The main idea of this approach 



can be explained as follows. Tunneling can be represented as a result of the accumulation 
of quantum fluctuations whose amplitude greatly exceeds their usual value determined by 
the uncertainty principle. This happens when the long wavelength quantum fluctuations 
responsible for the tunneling correspond to bosonic excitations with large occupation num- 
bers. In such cases one can treat these fluctuations as classical fields experiencing Brownian 
motion due to their interaction with the short wavelength quantum fluctuations. 

Suppose that the large fluctuations of the scalar field responsible for reheating in the 
model (|27|) initially look like spherically symmetric bubbles (which is the case if the proba- 
bility of such fluctuations is strongly suppressed, see above). The equation of motion for a 
bubble of a scalar field 0(r) in Minkowski space is 

= 0" + 20V-1 - r'(0) . (30) 

Here r is a distance from the center of the bubble and 0' = At the moment of its 
formation, the bubble wall does not move, = 0, = (critical bubble). Then it gradually 
starts growing, > 0, which requires that 

|0" + 20V-i| < -K'(0) . (31) 

A bubble of a classical field is formed only if it contains a sufficiently large field 0, and if the 
bubble itself is sufficiently large. If the size of the bubble is too small, the gradient terms 
are greater than the term |y(0)|, and the field inside the bubble does not grow. 

At small r the shape of the bubble can be approximated by = 0(0) — «r^/2. In this 



approximation, the bubble has a typical size tq ~ y ^^7^? and 0'r ^ = 0" = —a. Therefore 
at the moment of the bubble formation, when = 0, one has 
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0" = r(0(O))/3 



(32) 



Replacing 0" by A;q0(O) one finds that the bubble can be considered a result of overlapping 
of quantum fluctuations with typical momenta k ^ ko r^^ , where 



,.^C^-nm), (33) 

° 30(0) ^ ^ 

Here C = 0(1) is some numerical factor reflecting uncertainty in our estimate of fco- 

Let us estimate the probability of an event when vacuum fluctuations occasionally build 
up a configuration of the field satisfying this condition. In order to do it one should remember 
that the dispersion of quantum fluctuations of the field with k < ko'is given by (50^) ~ 
This gives 

,.,2. ^2 V'im) .... 

{^^ )^<^o -^2=C . (34) 

This is an estimate of the dispersion of perturbations that may sum up to produce 
a bubble of the field that satisfies the condition (|3T|). Of course, this estimate is rather 
crude. But let us nevertheless use eq. (0) to evaluate the probability that these fluctuations 
build up a bubble of a radius r ^ k^^ containing the field at its center. Assuming, in the 
first approximation, that the probability distribution is gaussian, one finds: 

Let us first apply this result to the theory — A0'^/4. In this case one finds 

n0)~exp(^-^) . (36) 

Note that the factor in the exponent in (|36| ) to within a factor of C ~ 2 coincides with the 
Euclidean action Se in eq. (|^). Taking into account the very rough method we used to 
estimate fco and calculate the dispersion of the perturbations responsible for tunneling, the 
coincidence is rather impressive. It was shown in p2| , p!9| that this approach gives exactly 
the same answer as the Euclidean approach for the case of tunneling during inflation when 
V" < H^. 

Most importantly, this method allows us to investigate tunneling and the development 



of instability in the theories where instanton solutions do not exist |T9|. In particular, for 
tunneling in the theory — Af 0'^/3 one finds 

P(0)~(A.;0)^exp(^-^) . (37) 

We included here the subexponential factor 0(A;q) ~ (Af0)^, which is necessary to describe 
the probability of tunneling per unit time per unit volume. 

This means that tunneling is not suppressed for ~ This result is in agreement 

with our previous estimate (pSf). Now let us take into account that the total time of the 
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development of instability is a sum of the time of tunneling plus the time necessary for rolling 
of the field down. One can show that the time of rolling down is inversely proportional to 
m(0) ~ i/Af 0, i.e. it decreases at large 0. Also, the subexponential factor {XvcpY grows 
at large 0, which makes tunneling to large faster. Consequently, as we already discussed 
above, the main contribution to the development of instability is given by the fluctuations 
with ^ j2^- Exponential suppression of the probability of such fluctuations leads to 
their approximate spherical symmetry. 

The results of our lattice simulations for this model are shown in Fig. |12|. In this model 
bubbles form quickly enough (unlike in the model — A0^), so we were able to start with 
quantum fluctuations centered at (0) = and allow the bubbles to form. The bubbles (high 
peaks of the field distribution) grow, change shape, and interact with each other, rapidly 
dissipating the vacuum energy V{0). 
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FIG. 12. Field values on a 2D slice through the lattice for V = -^v(t)^ + f c/)"^ The 



growth of quantum fluctuations of (p looks like bubble formation. Remarkably, the bubbles expand 
and collide even before the average field value reaches the minimum. Preheating occurs due to a 
combined effect of bubble production, tachyonic instability and bubble wall collisions. This figure 



should be compared with Fig. |^ for the theory V 



+ 1 . 
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Figure |T3| shows the probabihty distribution P(0, t) in the model (]27|). As we see, in this 
model the field distribution also rapidly relaxes near the minimum of the effective potential 
within a single oscillation. In this case the histogram in the beginning looks pretty chaotic 
because of bubble formation and bubble wall collisions, which we could see in the previous 
figure. 
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FIG. 13. Histograms describing the process of symmetry breaking in the model ( p7[ ) for 
A = 10^^. After reaching the minimum of the effective potential, the distribution acquires the 
form shown in the last frame and practically does not oscillate. The last histogram corresponds to 
the last frame in Fig. 

One should note that the numerical investigation of this model involved specific compli- 
cations due to the necessity of performing renormalization. Lattice simulations involve the 
study of modes with large momenta that are limited by the inverse lattice spacing. These 
modes give an additional contribution to the effective parameters of the model. In the limit 
of zero lattice spacing these corrections would become infinite, but they are regularized by 
the lattice cutoff. In our simulations of the simple model ([^) these corrections gave a contri- 
bution to the effective mass of the field that was much smaller than m and therefore did 
not affect our results. Meanwhile, in the cubic model similar corrections induce a (fictitious) 
linear term Af0(0^). This term should be subtracted by the proper renormalization proce- 
dure, which brings the effective potential back to its form (P7|). See the appendix for more 
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details. This was the first time in our simulations when a careful treatment of high frequency 
modes was necessary. A similar situation may occur in any theory where V"{0) = 0, such 
as the theory —X(j)'^ discussed in the previous section. 

Figure |l^ shows the occupation numbers of produced particles in the model (^) with 
A = 10^". These occupation numbers grow up to 10^ - 10^ within a single oscillation, which 
is in good agreement with our estimate (pOl). 




IX. CONCLUSIONS 

In this paper we studied the dynamics of spontaneous symmetry breaking, which occurs 
when a scalar field falls down from the top of its effective potential. We have found, in 
agreement with [|^, that the main part of this process typically completes within a single 
oscillation of the distribution of the scalar field. This is a very unexpected conclusion that 
may have important cosmological implications. 

One of the most efficient mechanisms for the creation of matter after inflation in theories 
with convex effective potentials {V"{(j)) > 0) is the mechanism of parametric amplification 
of vacuum fluctuations in the process of homogeneous oscillations of the infiaton field, which 
was called preheating M. It has also been noted that in the case where potentials become 



concave {V"{(f)) < 0), preheating may become more efficient |31]. Now we see that this effect 
is very generic. In many theories with concave potentials the energy of an unstable vacuum 
state is transferred to the energy of inhomogeneous classical waves of scalar fields within a 
single oscillation of the field distribution. We emphasize here that we are talking about the 
oscillations of the field distribution rather than about the oscillations of a homogeneous field 
because quite often the homogeneous component (0) of the field remains zero during 
the process of spontaneous symmetry breaking. 
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One of the important consequences of our results is the observation [|l| that in many 
models of hybrid inflation |3^ the first stage of reheating occurs not due to homogeneous 
oscillations of the scalar field but due to tachyonic preheating A detailed discussion of 
this effect will be contained in 0. 

The process of preheating and symmetry breaking may take an especially unusual form in 
the theory of brane inflation |3^j3^ based on the hybrid inflation scenario and the mechanism 
of tachyon condensation on the brane antibrane system . 

The situation in models of the type used in the new inflation scenario is somewhat more 
complicated. In these models the potential is also concave. However, the expansion of the 
universe stretches inhomogeneities of the field rolling down from the top of the effective 
potential and makes it homogeneous on an exponentially large scale. Therefore to evaluate 
a possible significance of tachyonic instability in this regime one must compare the amplitude 
of the homogeneous component of the field with the amplitude of the quantum fluctuations. 
The result appears to be very sensitive to the scale of spontaneous symmetry breaking in such 
models. A preliminary investigation of this issue indicates that in small-field models where 
the scale of spontaneous symmetry breaking is much smaller than Mp, the leading mechanism 
of preheating typically is tachyonic. If correct, this would be a very interesting conclusion 
indicating that in large-field models the leading mechanism of preheating typically is related 
to parametric resonance, whereas in small-field models the main mechanism of preheating is 
typically tachyonic, at least at the first stages of the process. We will return to the discussion 
of this issue in a coming publication p. 

Finally we should mention that an interesting application of our methods can be found in 
the recently proposed ekpyrotic and pyrotechnic scenario [^,^]. Even though we are very 
skeptical with respect to the ekpyrotic/pyrotechnic scenario for many reasons explained in 
P7| , it is still interesting that the methods developed in the theory of tachyonic preheating 
provide us with a very simple theory of the generation of density perturbations in these 
models [P7[]. 

It is a pleasure to express our gratitude to J. Garcia-Bellido, P. B. Greene, and I. Tkachev 
for the collaboration at the early stages of this project. We thank NATO Linkage Grant 
975389 for support. L.K. was supported by NSERC and by CIAR; G.F. and A.L. were 
supported by NSF grant PHY-9870115 and by the Templeton Foundation. 



APPENDIX A: THE LATTICE CALCULATIONS 
1. Overview 

The lattice calculations reported on in this article were all done using the program 
LATTICEEASY, developed by Gary Felder and Igor Tkachev. The program records the 
value of the fields and derivatives at each point on a spatial grid with evenly spaced points. 
The fields are then evolved using their classical equations of motion 

- V|0 +V' = 0. (Al) 

The use of the classical equations is justified because the instability discussed rapidly drives 
the fields to a state with exponentially large occupation numbers, meaning they effectively 
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act as classical fields |jTO|]. Although LATTICEEASY is designed to (optionally) include 
the effects of cosmological expansion on field evolution all of the simulations reported here 
were done in a fiat spacetime background. The effects of expansion will be discussed in 
subsequent publications ( @,|]). 

Time evolution is done with a staggered leapfrog algorithm using a fixed time step. The 
initial conditions for the fields and derivatives are set in momentum space and then Fourier 
transformed to give the initial spatial distribution. The initial values of the modes are given 
by quantum fluctuations. Each mode has a random phase and a Gaussian random amplitude 
with expectation value 

m') = ^- (A2) 

The exception to this is the Coleman- Weinberg potential for which the initial conditions 
were set by the instanton configuration described in the text. Note that ordinarily the 
equations for quantum fluctuations would have yW^+rri? in the place of k in the formula 
above. For the models discussed here, however, we were simulating a quench in which the 
effective squared mass of the fields is presumed to have rapidly become negative. Thus we 
used initial conditions corresponding to massless fluctuations. For some of the runs here we 
imposed a momentum cutoff, setting 0^ = for all modes above a certain momentum k. 
Such cutoffs eliminated unphysical effects from quantum fluctuations that were not excited 
to large values. In each such case we also ran without the cutoff and found the results to be 
qualitatively similar except for the addition of high frequency noise in the field distribution. 

The plots shown in this paper show either field values (which are self-explanatory), 
probability distribution functions (PDF), or occupation number spectra. The PDF of a 
field is obtained by dividing the field values on the grid into evenly spaced bins and simply 
counting the number of gridpoints in which the field value was in each bin. The occupation 
number is defined by Fourier transforming the field and computing for each mode by Eq. 
(0), where rn? = V", = ^/k"^ + m? for > 0. For < one can use either cuk = \k\ or 



LUk = yk"^ + |m^|. All plots shown in this paper except Fig. 1 and Fig. 2 use Uk = \k\ when 
< but the results were qualitatively similar using either definition. The occupation 

number is given by averaging over a spherical shell in Fourier space. This definition 

coincides with the standard one in the end of the process, where the mass squared becomes 

positive and topological defects disappear. 

The full details of these lattice calculations can be found in the documentation available 

on the LATTICEEASY website at |http:/ /physics. Stanford. edu/latticeeasy| . Moreover, these 



calculations have been discussed in previous publications of ours (e.g. ||TO|JlT[| ). This is our 
first publication where we discuss simulations that used renormalization, however, so we 
will discuss this procedure in the next section. The last section of the appendix lists the 
parameters used for each of the runs illustrated in the paper. 



2. Renormalization 

As we have discussed, the justification for doing a classical calculation for quantum fields 
is that once the field fluctuations are amplified sufficiently quantum effects are negligible. 
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There are some cases, however, when these quantum effects may be important, and in such 
cases they may be (partially) accounted for through a simple form of renormalization. 

Consider how this applies to the lattice calculations discussed here. Initially the field 
fluctuations are only those representing quantum vacuum states. These fluctuations affect 
couplings, masses, and the total energy of the system in a way that is dependent on the 
lattice spacing. For example, consider the theory 

V = ^-X<p' - \\v<^'' + ^^\v' (A3) 

and rewrite the field 0(a;, t) as the sum of a homogeneous component (/)(t) and fluctuations 
50. The effective potential felt by the homogeneous field will receive a correction (among 
others) from the fluctuations equal to 

bV ^ -Xv{6(j)'^)(f). (A4) 

(The 1/3 is cancelled by a coefficient arising from combinatorics.) This correction represents 
an unphysical effect in the sense that its strength depends on the ultraviolet cutoff imposed 
by the lattice. In the limit of zero lattice spacing where arbitrarily large momenta would be 
included on the lattice this correction would become infinite. This would add an unphysical 
term C0 to the effective potential. This effect can be eliminated, however, by adding a 
counterterm 

AV = Xv{6(j)'^)(f) (A5) 

or equivalently by adding the term 

Av(50') (A6) 

to the equation of motion for 0. Note that (50^) in this case refers to the value that arises 
from initial quantum fluctuations, not to a dynamic quantity that changes as the field evolves 
and fluctuations grow. Such changes represent physical effects and should not be eliminated. 
In effect this correction eliminates the linear term in the potential at = when the field 
is in the vacuum state. 

The above example illustrates how a simple form of renormalization can be implemented 
on the lattice. This procedure could in principle be used to renormalize any mass, coupling 
constant, or energy term in the theory. Ordinarily these corrections are not important 
(unless one uses a very large value of the momentum cutoff) because the quantum effects 
are quickly swamped as the fluctuations become amplified. We did not find it necessary to 
use renormalization for any models except the cubic one, where we used it as described here 
to prevent the field from artificially rolling away from = due to the induced linear term. 

3. List of Parameters 

In this section we list the parameters used for the lattice simulations from which all of 
the figures in the paper were drawn. The models discussed in the paper will be referred to 
here simply as Quadratic (model 0), Complex (model ^ with 0^ replaced by |0P), Quartic 
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(model 0), and Cubic (model |2^). The parameters A^, L, dt, and kcut refer to the number of 
gridpoints, the width of the box, the time step, and the initial momentum cutoff respectively. 
Lengths and times are measured in units of y/Xv. The coupling constant A is also given for 
each run, as well as any information specific to the particular plot. All runs are assumed to 
be three dimensional unless otherwise indicated. 



Figures |, |, and | Quadratic model, = 256^, L = 100, dt = .1, k^ut = 1, A = lO"! 
The run shown in Fig. ^ started with a homogeneous displacement (0) = .If. The 
PDF's in figure ^ and |^ use 256 bins. 

Figure |^ Quadratic model. Two dimensional simulation with = 1024^, L = 1600, 
dt = .1, kcut = 0, A = 10~^. The slices shown are plots of 128^ points where each 
point represents an average over an 8 x 8 box of gridpoints. This averaging was done 
simply to make plotting easier and reduce the resulting image sizes. 

Figures ^ and ^ Complex model, = 256'^, L = 600, dt = .2, A = 10^^. The momentum 
cutoff was kcut = 2 and kcut = for Figs. |^ and ^ respectively. The PDF's in Fig. |^ 
used 80 X 80 bins for the 2D field space. The spectra shown in Fig. ^ are for the real 
component Re(0). The spectra for lm(0) look identical. 
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Figure |^ Complex model, A^ = 128^, L = 40, dt = .01, kcut = v^, A 

defined as the collection of points on the lattice for which |0p < .02. 



-4 



A string was 



Figures 10 and O Quartic model, A^ = 256^ L = 500, dt = .1, 



Both runs used as initial conditions the instanton solution (P 



kcut 
at t 



0, A = 
with 



5 X 10-5. 
^0 = .02v. 



The run shown in Fig. ITOl did not include additional quantum fluctuations, while the 



run shown in Fig. |TT] did. Both figures show a partial 2D slice through the lattice, cut 
so as to show the interior of the bubble. 



10" 



The 



Figures [T|, and |T4| Cubic model, A^ = 256^ L = 200, dt = .25, A 

momentum cutoff kcut was .6 for Figs. |12| and |T3| and for figure 0. The PDF*s in 
Fig. 



13 used 256 bins. 
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